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\^ , We report the results of a blind search for 22-GHz water masers in two regions, 

fSJ ' covering approximately half a square degree, within the giant molecular cloud asso- 

I dated with RCW 106. The complete search of the two regions was carried out with 

■ the 26-m Mount Pleasant radio telescope and resulted in the detection of nine water 

' masers, five of which are new detections. Australia Telescope Compact Array (ATCA) 

r~| . observations of these detections have allowed us to obtain positions with arcsecond 

^ i' accuracy, allowing meaningful comparison with infrared and molecular data for the 

region. We find that for the regions surveyed there are more water masers than ei- 
ther 6.7-GHz methanol, or main-line OH masers. The water masers are concentrated 
towards the central axis of the star formation region, in contrast to the 6.7-GHz 
^ _ methanol masers which tend to be located near the periphery. The colours of the 

GLIMPSE point sources associated with the water masers are similar to those of 6.7- 
GHz methanol masers, but slightly less red. We have made a statistical investigation 
of the properties of the ^^CO and 1.2-mm dust clumps with and without associated 
^ ■ water masers. We find that the water masers are associated with the more massive, 

' denser and brighter ^'^CO and 1.2-mm dust clumps. We present statistical models 

that are able to predict those ^'^CO and 1.2-mm dust clumps that are likely to have 
associated water masers, with a low misclassification rate. 
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1 INTRODUCTION 

Water masers have been regarded as one of the best indica- 
tors of star formation since soon after the first detection of 
22-GHz water mase r emission in the sou rces Sgr B2, Orion 
and W 49 in 1969 feheun g et allll969l ). The 22-GHz 6i,6 
52,3 rotational transition of H2O is the brightest spec- 
tral line at radio wavelengths and traces shocked gas in star 
formation regions, outflows, as well as dense circumstellar 
shells around evolved stars. Emission from this transition of- 
ten exhibits significantly greater temporal variability than is 
commonly observed i n other interstellar maser species such 
as OH and methanol l|Brand et aLlbOO^ ). The physical con- 
ditions required to produce water maser excitation are high 
densities (10^-10^ cm~^) and tempe ratures of a few 100 K 
lElitzur. HoUenbach fc McKeelll989l ') , both of which are seen 



in the inner parts of circumstellar disk s around young stella r 
objects and in regions of shocked gas l|Torrelles et al.ir2002l ). 

Interstellar ma sers from wate r (e.g. IValdettaro et al] 

200l[). OH (e.g. ICaswelll 119981 ) and methanol (e.g. 



Pestalozzi. Minier fc BoothI 20051 ) transitions have been de- 



tected towards hundreds of star formation regions in our 
Galaxy, with many of these showing emission from more 
than one species. While there have been a number of large- 
scale untargeted sear ches for OH (ICaswell. Havnes fc Goss 
1980) and methanol (|Ellingsen et al.lll996l : ISzvmczak et al 
20021) ■ previous searches for water masers in star forma- 



tion regions have typically targeted ultracompact Hll regions 
selected on the b asis of IRAS (Infrared Astronomy Sate l- 
lite) colours fe.g. IChurchwell. Walmslev fc Cesaronilll990l ). 
or other s ources believed to be high-mass y oung stellar ob- 
jects (e.g. iBeuther. Walsh fc SchiM |2002| ) . To date there 
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have been no large untargeted surveys for water masers pri- 
marily because at a frequency of 22 GHz telescope beam 
sizes are approximately one-third the size of those at 6.7- 
GHz, the frequency of the strongest methanol maser tran- 
sition, and hence require approximately an order of magni- 
tude more pointings. Here we present an untargeted search 
for water masers within the giant molecular cloud (CMC) 
complex associated with RCW 106 (the G 333.2-0.6 GMC). 
The GMC is located a t a distance of 3 . 6 kpc (|Lockman 19791 ) 
and was discovered bv I Gillespie et ahl l|l977h during observa- 
tions of molecular clouds associated with southern Galactic 
Hll regions in the J=l-0 transition of CO. These observa- 
tions uncovered a number of bright Hll regions along a line 
which is almost parallel to the Galactic plane including one 
of the brightest infrared s ources in our Galaxy, IRAS1Q1%2>— 
4958 l|Becklin et al.lll973l ). which is associated with the Hll 
region G 333.6-0.2. 

The GMC is roughly centred on / ~ 333°, 6 ~ -0?5 (or 
Oj20oo = 16:21, (5j2ooo=-50:20) and extends approximately 
1?2 X 0?6 on the sky (or approxi mately 90 pc x 3 pc at 
an assumed distance of 3.6 kpc ([Bains et aLll2006l )). This 
region passes through the ring of molecular clouds that 
circle the Galaxy at around 5 kpc from its centre (e.g. 
ISimon et al.l 1200 ll ) and exhibits a diverse range of molec- 
ular regions, bright Hll regions, GLIMPSE {Spitzer Galac- 
tic Legacy Infrared Mid-Plane Survey Extraordinaire) point 
sources, IRAS and MSX (Midcourse Space Experiment) 
sources, all of which are embedded in a larger region of dif- 
fuse atomic and molecular gas. 

Obs ervations of the ^ ^CO J=l-0 transition at 110 
GHz by I Bains et al.l (|2006l ) showed the velocity structure 
of the region to contain five primary velocity components, 
with the dominen t feature centred on ulsr ~ -50 km s~^. 
iBains et al.l (|2006l l showed that it is likely that at least five 
distinct molecular clouds lie along the line-of-sight (see Sec- 
tion [3XTJ. For this reason the G 333.2-0.6 GMC is ref- 
ered to as the main cloud, showing emission over the ve- 
locity range -35 km s~^ to -65 km s~^. Analysis of the 
integrated ^^CO d ata using the CLUMPFIND algorithm of 
IWilliams. de Geus fc Blitz (1994) identified 61 ^^CO clumps 
within the main cloud. The ^^CO emission takes the form 
of a string of knots with the clumps a rranged along an axis 
aligned NW to SE (|Bains et al.ll2006l ). 

This GMC has been the focus of numerous observa- 
tions in recent times, including far-infrared (FIR) obser- 
vations of the dust continuum at 150 and 210 /im which 
identified 23 emissio n peaks with dust temperatures be- 
tween 20 and 4 K jKarnik et al" 120011 ). The region was 
also observed bv lMookeriea et a" 1 20041 ) who used SIMBA 
(SEST IMaging Bolometer Array) on the Swedish European 
Southern Observatory Submillimetre Telescope (SEST) to 
obtained a 1.2-m m dust continuum image of the region. 
iMookeriea et al" I {2004) identified 95 dust emission peaks 
(or clumps), half of which have MSX counterparts. Obser- 
vations of a multitude of molecular lines towards detected 
Hll regions and IRAS sources indi cate probable ongoing star 
formation (Mookcriea et al.l l2004l ). These observations, like 
those of Bains ct a,l. (2 006. ) have identified that the GMC 
has a linear clumpy structure. 

Comple te surveys of th e GMC region have been car- 
ried out bv ^EUingsen et al] (Il996|) for 6.7 -GHz methanol 
masers and ICaswell. Havnes fc Gos3 
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Figure 1. Integrated ^"^CO emission observed by I Bains et al.l 
1I2OO6I ) with the two regions surveyed for water masers overlaid. 



1667-MHz OH masers. These surveys resulted in the detec- 
tion of nine methanol and six OH ma s ers w ithin the re- 
gion surveyed in '^^CO bv iBains et al.l (|2006l ). Five water 
masers have also bee n detected withi n this r egi on in targeted 
search e s made by Johnston et al.l (1 19721). Caswell et al.l 
(ll974h . lBatchelor et al.l (jl980') and'Braz & Scalisd (|l982l ). 

Here, we present the results of an untargeted survey of 
22- GHz water masers made with the University of Tasma- 
nia Mount Pleasant 26-m radio telescope. The survey cov- 
ers two distinct regions within the GMC. The first (here- 
after Region 1) is a 0?50 x 0?42 area centered on I ~ 333, 
b ~ -0?5. This region encompasses m uch of the hi g h den - 
sity gas and dust regio ns ide ntified by Karnik et al.l (|200lf ) , 
IMookeriea et al.1 (|2004h and iBains et al.1 (|2006l ) in the cen- 
tral section of the GMC. Region 1 contains one previously 
detected 22-GHz wa ter maser, G 333.13-0.43, discovered by 
ICaswell et all \l974 ). The second region (hereafter Region 
2) covers a 0?28 x 0?24 area of the GMC and is approx- 
imately centred on well-known optically visible Hll region 
RCW 106. This region contains two previously detected wa- 
ter masers. The extent of the two regions compared to the 
integrated ^^CO emission are shown in Fig. [1] 

The GMC is the focus of an ongoing project to char- 
acterise the turbulence in the molecular cloud and compare 
this to the star formation efficiency in order to attain a re- 
lationship between the two. Commencing in 2004 a multi- 
tude of millimetre molecular line transitions (including ^"^CO 
jBains et al.ll2006[ ). CS, C^^S, C^^O, C2H, HON, H^^^CN, 
HCO+, H^='CO+, HNC, CH3OH and SiO) have been ob- 
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served by the Delta Quadrant Survey Team at the University 
of New South Wale^ (UNSW). Interstellar masers require 
special physical conditions and the different species are gen- 
erally thought to trace particular evolutionary phases of the 
high-mass star formation process. Through combining in- 
formation on all the strong and common interstellar maser 
transitions with millimetre molecular line and other exist- 
ing infrared and millimetre continuum data we aim to better 
understand the evolutionary phases traced by each type of 
maser. 



2 OBSERVATIONS & DATA PROCESSING 

The primary search for 22-GHz water masers towards the 
G 333.2-0.6 GMC was undertaken between 2005 April- 
November using the University of Tasmania's 26-m radio 
telescope at Mount Pleasant. The observations were made 
with a cryogenically cooled receiver that detects both left 
and right circularly polarized signals and has a typical sys- 
tem equivalent flux density of 2000-2200 Jy in good weather 
at elevations above about 40 degrees. At 22 GHz the tele- 
scope has a 2.2-arcmin half power beam width (HPBW) and 
at the time of the observations the measured RMS pointing 
errors were ~1 arcminute. The data were recorded using 
a 2-bit auto-correlation spectrometer configured with 2048 
spectral channels per polarisation over a 32-MHz bandwidth, 
which was centred on an LSR velocity of approximately 
-40 km s~^. The observations covered a velocity range of 
430 km s^^ with a spectral resolution of 0.25 km s~^. The 
two regions were surveyed using an equilateral triangle pat- 
tern with each pointing separated by 1.1 arcminutes (half 
of the HPBW) from all adjacent grid points. Each point- 
ing was observed for a total onsource integration time of 
10 minutes. The water maser G 333.608-0.215, discovered 
bv I Johnston et af] l| 19721) . was observed at the beginning of 
each observing session to test the system and ensure con- 
sistency. Although it is not located in either of the survey 
regions we have included it in the nine masers detected in 
this survey. 

The weather conditions in which the observations were 
made varied substantially, and in general, the data for Re- 
gion 1 were taken under much more favourable conditions. 
In order to minimize atmospheric and telescope pointing ef- 
fects all observations were made above an elevation of 40 
degrees. Data affected by poor weather conditions were re- 
observed (in some cases on multiple occasions), however a 
small percentage of the data has a significantly poorer sensi- 
tivity. The observations of Region 1 required approximately 
750 pointings while Region 2 was much smaller containing 
just under 300 pointings. After box car smoothing over 5 
channels and averaging the two polarisations together the 
majority (78 percent) of the data taken in Region 1 had an 
RMS noise level of less than 1 Jy while the remaining (22 
percent of) data in this region was subject to an RMS noise 
level of between 1 Jy and 2 Jy. This equates to a 5(j detection 
limit of 5 Jy over the majority of the region, with a maxi- 
mum of 10 Jy in the worst affected pointings. In Region 2 
only 35 percent of pointings have an RMS noise level of less 

^ http:/ /www. phys.unsw.edu.au/astro/mopra/dqs.php 



than 1 Jy and the remaining 65 percent had an RMS noise 
of up to 5 Jy. This equates to a 5(t detection limit between 5 
and 25 Jy, hence our ability to detect weaker masers within 
Region 2 is greatly reduced. 

For each detected water maser additional observations 
consisting of a 5-point grid (centred on the preliminary posi- 
tion) were made to locate the emission more accurately. The 
position was determined by fitting a 2D circular Gaussian 
(with the same HPBW as the telescope) to the relative am- 
plitudes of the strong emission in the maser spectra observed 
in the 5-point grid. Positions determined in this manner are 
accurate to approximately 1 arcminute, which is insufficient 
to allow meaningful comparison with millimetre and infrared 
observations of the region. In order to obtain accurate posi- 
tions for the detected water masers we were granted two ses- 
sions of Australia Telescope Compact Array (ATCA) direc- 
tor's time. Preliminary observations of the masers detected 
in the Mount Pleasant survey were made on 2006 June 19 
with the ATCA in the 1.5D configuration. The observations 
were centred on a frequency of 22.238 GHz with the cor- 
relator configured to record 256 spectral channels across a 
16-MHz bandwidth. These observations failed to detect one 
of the water masers discovered in the Mount Pleasant survey 
(G 333.29-0.38), most likely due to temporal variability. 

Further observations were made with the ATCA in the 
6A configuration on 2006 July 16 & 17. In this array config- 
uration the minimum baseline length is 337 m and the max- 
imum is 5939 m. The observations were centred on 22.236 
GHz and the correlator sampled two orthogonal linear po- 
larisations, each processed to give a 512 channel spectrum 
across an 8-MHz bandwidth. Each of the maser sites de- 
tected in the earlier ATCA observations (8 in total) were 
observed in a series of three minute cuts over a range of 
hour-angles. Observations of one of the phase calibrators 
MRC B1613-586 or MRC B1646-50 were made for a dura- 
tion of 90 seconds before and after every two maser observa- 
tions (i.e. every 6 minutes). PKSB1934-638 was used as the 
primary flux calibrator and at 22.236 GHz has an assumed 
flux density of 0.83 Jy. PKSB1253-055 was used for band- 
pass calibration. These ATCA observations were made using 
reference pointing, and as arcsecond accurate positions were 
obtained in the initial ATCA observations all of the maser 
sites were close to the centre of the primary beam. Over the 
two days each source was observed a total of eight times, 
equivalent to a total on-source integration time of 24 min- 
utes. 

The data were processed using t he miriad software 
package (jSault. Teuben fc Wrightlll995r ) applying the stan- 
dard techniques for ATCA spectral line and continuum ob- 
servations. Continuum sources which were located away 
from the centre of the primary beam have had their flux 
densities corrected to account for beam attenuation. The 
frequency resolution, after Banning smoothing, was 0.038 
MHz or 0.50 km s'^ The RMS noise in a single spectral 
channel in the final data cubes was approximately 0.15 Jy 
and the signal-to-noise ratio of the final spectra was about 
20:1 in the worst case. The RMS noise level in the continuum 
images was typically around 0.02 Jy/beam. These observa- 
tions have enabled us to determine the positions of the water 
masers to an accuracy of approximately 0.5 arcseconds. 
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3 RESULTS 

A search of two regions near the G 333.2-0.6 giant molec- 
ular cloud resulted in the detection of nine 22-GHz water 
masers, five of which are new detections (Table as well 
four 22-GHz continuum sources (see Section [3.2p . Figure [2] 
shows the positions of all the detected water ma sers overlay- 
ing th e integrated ^^CO emission observed by iBains et al.l 
l|2006l ). while Fig. |3] shows the maser locations on a three 
colour GLIMPSE image of the GMC. Comparison of the 
water maser locations with the other maser species, the inte- 
grated ^^CO emission and the three colour GLIMPSE image 
shows that in general the water masers originate very close 
to the higher density molecular gas and warm dust, near 
the main axis of star formation within the molecular cloud. 
In contrast the methanol masers tend to be offset from this 
axis, close to the interface between the i ntense mid-infra red 
emission and the larger molecular cloud (|Ellingsenll2006l l . 

Spectra of the detected water masers are shown in 
Figs. U and [S] The spectra have been produced by integrating 
the emission in the ATCA image cubes for each source. The 
only exception is water maser G 333.29-0.38 for which the 
Mount Pleasant spectrum is shown. The positional accuracy 
of sources detected in the ATCA observations is approxi- 
mately 0.5 arcseconds and we have used three significant 
figures after the decimal place in their Galactic coordinate 
names. For the source only detected in the Mount Pleasant 
component of the survey we are only justified in using 2 
significant figures and have done so throughout the paper. 
Comments on each maser can be found in section 13.11 The 
newly detected water masers (Fig. |4]), with one exception 
have a peak flux density less than 50 Jy, while the previ- 
ously detected sources (Fig. again with one exception, 
have peak flux densities greater than 100 Jy. 

The 6.7-GHz methanol masers sites in this region have 
previously be en search e d for a ssociated 22-GHz water maser 
emission by iHanslowl who detected emission to- 

wards a number of sources (G 332.942-0.686, G 333. 121- 
0.434, G 333.128-0.440, G 333.130-0.560, G 333.234-0.062 
and G 333.466-0.164) in the G 333.2-0.6 giant molecular 
cloud. Three of these (G 332.942-0.686, G 333.234-0.062 and 
G 333.466-0.164) lie outside the regions of our untargeted 
search. While of the other three, only G 333.121-0.434 was 
detected in the current work. It appears that the emis- 
sion attributed to G 333.128-0.440 by Hanslo^ (119971 ) is in 
fact associated with G 333.121-0.434, while that associated 
with G 333.130-0.560 was not detectable at the epoch of our 
search. Because of the uncertaint y in th e position of the wa- 
ter masers detected by IHanslowl |l993) and the possibility 
that some may result from an unassociated strong source 
detected in a sidelobe, in Figs. [2] & 13 we have only marked 
those water masers detected in our current survey. 

In addition to the nine water masers that we present, we 
made a one-time detection with the 26-m Mount Pleasant 
radio telescope of a water maser right on the edge of the field 
of observations with coordinates G 333.22-0.20. Subsequent 
observations showed no detectable emission and as a result 
we do not include it in the list of water masers that we 
detect. This emission consisted of a single velocity feature at 
-87 km s~^ of around 6 Jy. We believe that this emission was 
actually a detecti on of the stron g water maser G 333.234- 
0.062 detected bv IHanslowl l| 19971 '). IHanslowl l| 19971 ') reported 



G 333.234-0.062 to consists of multiple velocity features with 
the most prominent observed at -86 km s~^, with a fiux 
density of 108 Jy. 

The majority of the water masers we detected have ex- 
hibited variability of up to a factor of 10 on a time scale of 
several months. This type of variability is common in water 
masers and a survey of water maser emissio n towards main- 
line O H masers in star formation regions bv lBatchelor et al] 
l|l980l ) found that about 60 percent of water masers exhib- 
ited variability of up to a factor of two over an eight month 
period, while the remaining 40 percent exhibited more ex- 
tensive variability. Given that our initial observations were 
made in varying weather conditions with comparitively poor 
pointing accuracy it is difficult to quantify the absolute vari- 
ations accurately. However, because the water masers have 
multiple spectral components we are able to determine that 
variability has occurred by examining the relative ampli- 
tudes. Water masers associated with low mass stars are typ- 
ically both weaker and more variable than t hose associated 
with high mass stars (|Claussen et al.lll996l ). So we would 
expect our observations to be more likely to detect masers 
associated with high-mass star formation, than those asso- 
ciated with less luminous objects. Variability is also a likely 
explanation as to why two of the masers detected in the 
survey with the 26-m Mount Pleasant radio telescope were 
not detected in the final observations made with the ATCA 
6-8 months later, particularly if these masers are associated 
with low-mass stars. 

We have compared the positions of the water masers 
that we detected with GLIMP SE, IRAS and MSX so urces, 
as well as 1.2- mm dust clumps l|Mookeriea et all 20041'). FIR 
sources (|Karnik et al.1 l200ll ). " CO emission (iBains et al.l 
l2006h . CS emission (N. Lo, private communication) and 
other maser species. The relative positional accuracy of each 
of these datasets differs, some having significantly poorer po- 
sitional accuracy than our ATCA observations. We consider 
a water maser to be associated with a GLIMPSE, IRAS, 
MSX or FIR source if it is within a radius of 2, 30, 5 or 60 arc- 
seconds respectively. For water maser G 333.29-0.38 we use 
more relaxed positional constraints as its position is much 
less accurately known. For this source we consider the water 
maser to be associated if it is within 5, 90, 40 or 60 arcsec- 
onds of a GLIMPSE, IRAS, MSX or FIR source respectively. 
We consider a water maser to be associated with a ^^CO, CS 
or 1.2-mm dust continuum clump if its position falls within 
the radius of the ^^CO, CS or 1.2 m m dust clum p s. We 
have use d the clump radii reported by iBains et al.l (|2006l ) 
and Moo keriea et al.l l\2004 ) for the ^^CO and 1.2-mm dust 
clumps, respectively. In the case of the CS data we have 
considered that a water maser falls within the radius of a 
given CS clump if it lies within emission that has a value of 
fiux density per beam which surpasses a 5a level. Velocity 
ranges of the ^^CO and CS clumps correlate well with the 
velocity ranges of the associated masers (see Tables [7] and 
[5]). Table [2] summarises all associations. 

Of the nine water masers that we detect, four are as- 
sociated with a GLIMPSE point source, three are associ- 
ated with an IRAS source, two are associated with an MSX 
source, five are associated with a FIR clump (|Karnik et al.l 
booiT) and seven are as sociated with a 1.2-mm dust clump 
( Mookeriea et al.ll2004 ). All of the water masers detected 
in this survey either lie within a ^^CO clump identified by 
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iBains et alj (|2006l ) or an identifiable emission peak of tfie 
^^CO data. In addition all of the water masers, for which 
CS data of the GMC was available, are associated with CS 
emission peaks covering a comparable velocity range to the 
masers. 

3.1 Individual sources 

Comments on the nine water masers that we detect can be 
found in Section 13.1.21 Refer to Table [1] for a summary of 
each of the sources. Figures [2] and |3] show the positions of 
the water masers overlaying the integrated ^''CO emission 
iBains et al.ll2006l ) and a three colour GLIMPSE image of 
the GMC (also showing the locations of 6.7-GHz methanol 
masers (jEllingson 2005.) and main-line OH masers ()CasweU 
Il998l )). respectively. 

3.1.1 Possible associations with a foreground cloud 

iBains et al.l (|2006l ) found that the mean velocity profile of 
the ^'^CO emission averaged over their field of observations 
exhibited five distinct velocity features (see Fig. [6]), while 
only the brightest velocity feature was attributed to the 
G 333.2-0.6 GMC. This most prominent feature is also the 
broadest feature and spans a velocity range of -35 km s~^ 
to -65 km s~^. The additional four velocity features seen in 
Fig. [6] a re centred on velo cities of -105, -90, -70 and -10 
km s" -\ iBains et al.l (|2006l ) plotted the kinematic distance 
versus the LSR velocity at the centre of the GMC and found 
that for the clouds at -10, -50, -70, -90, -105 km s"^ the 
associated distances were approximately 1, 3.5, 4.5, 5.5 and 
6.5 kpc. 

Water masers 1, 2, 4-7 and 9 have velocities which are 
comparab le to the broadest feature of velocity profile ob- 
served bv lBains et al.l l|2006l ) and as such we associate these 
masers with the main cloud. The velocity ranges of water 
masers 3 and 8 indicate that they are probably associated 
with a foreground cloud (or clouds) located at a distance of 
1 kpc. Water maser 3 is separated by less than an arcsecond 
from the 22-GHz continuum source G 333.060-0.489, indi- 
cating that the two are likely to be associated. This means 
that like water maser 3 this continuum source is likely to be 

associated with a foreground cloud. 

Visual inspection of the ^"^CO data of iBains et al.l 

l|2006D shows that at velocities of approximately -10 km s 
emission is restricted to a small range of right ascensions 
and declinations and is elongated almost perpendicular to 
the axis of the main GMC. Water maser 3 lies within this 
emission that is located at a distance of approximately 1 kpc. 
From the ^''CO data it appears unlikely that water maser 
8 is associated with the same cloud as the '^''CO emission 
associated with this maser does not overlap with the ^"^CO 
emission associated with water maser 3 and the spatial sep- 
aration between these masers is significant. 

3.1.2 Comments on individual sources 

1. G 332.653-0.621: Th is water maser was discovered by 
iKaufmann et al" 

1 l| 19761 ), who observed it to have a peak 
flux density of 58 Jy at a v elocity of -47 km s~^ in 1975. 
Subsequent observations by iBatchelor et al.l (|l980l ) during 



1977 May showed a single maser feature at -44 km s~^ 
and a slightly weaker intensity of 30 Jy. We measured the 
peak flux density to be 29 Jy at -45 km s~^, similar to the 
observations made in 1977. No OH or 6.7-GHz methanol 
masers have been detected associated with this water maser 
ICasweU. Havnes fc Gos3ll980l : lEllingsen et al.lll996l ). 

This so urce is offset from th e Hll region G 332.663-0.621 
identifled bv lHuang et al.l (|l999l ) by 35 arcseconds at a posi- 
tion angle of 58° and is situated within the RCW 106 com- 
plex. IRAS somce /T?^^ 16158-5055 is located 11 arcseconds 
away from this maser and has colours typical of an ultra- 
compact Hll region. This maser is separated from the centre 
of the nearest ^'^CO clump by 35 arcseconds and is 45 arc- 
seconds from the centre of the 1.2-m m dust emission peak 
identified by iMookeriea et al.) (|2004h at position angles of 
-79° and -83° respectively. 

2. G 332.826-0.5 49: This maser was discovered by 
iBraz fc Scahsd l|l982l) in 1980 April, who reported it to 
have a peak flux density of 250 Jy at -70.8 km s~^. We 
found the intensity peak of the maser to be 239 Jy at - 
59.1 km s~^ with emission covering a range of more than 
a 20 km s~^. While the velocity range over which emission 
has been observed has remained roughly constant since the 
maser's discovery, the relative intensities and velocity of the 
str ongest emission have not, for example observations made 
bv lBraz fc Scalisd l|l982l ) in 1981 May showed a peak at -62 
km s"^ of 198 Jy. 

This maser appears to be associated with the MSX 
source G 332.8269-00.5489 which is offset by 5 arcseconds 
and may be also associated with /T?^^ 16164-5046 which 
is located 29 arcseconds away, as well as a FIR source 
(|Karnik et al.|[2001iV The maser is probably a ssociated with 
a ^^CO clump identified bv lBains et all (|2006|) and a, 1. z-mm 
dust emission peak identified by Mookeriea et al.| l|2004l ). 
which are separated from the maser by 11 arcseconds and 
8 arcseconds respectively. This maser is located about 2 ar- 
cminutes from the peak of the RCW 106B complex which is 
centered on 16:20 -50:52 and is offset from the peak of the 
22-GHz continuum source detected in our ATCA observa- 
tions by 2 arcseconds. The w ater maser is s eparated from the 
OH maser G 332.826-0.548 l|Caswelllll993 ) by 8 arcseconds 
(see fourth sub- image of Fig. 

3. G 333.060-0.488: This maser exhibits several spectral 
features over an 18 km s~^ velocity range, with the most in- 
tense having a fiux density of 64 Jy at -8.7 km s"'^. This 
source was first observed at Mount Pleasant on 8 August 
2005 when the feature at -8.7 km s~^ was approximately 
6.5 Jy and the secondary feature at about km s~^ was 3 
Jy, implying a variation of a factor of 10 over an 11 month 
period. The velocity of this maser is comparable to one of 
the secondar y features of the velocity profile of the GMC 
observed by (|Bains et all 120061 ) (see Fig. [gj. This maser is 
likely part of a different molecular cloud located along the 
line of sight at a distance of approximately 1 kpc (see Sec- 
tion EXT)). 

The ATCA observations detected a 22-GHz contirmum 
source offset from the water maser source by about an arc- 
second (see third sub-image of Fig. [2]). The water maser 
source appears to be associated with GLIMPSE point source 
G 333.0600-00.4888. 

4. G 333 . 121-0 . 434: This source was discovered by 
ICaswell et al.l (|l974l ) in 1973 June. Subsequent observations 




Figure 2. The main image shows the integrated ^''CO emission observed by I Bains et ahl l|2006h with the positions of water (circle), 
methanol (square) and OH (cross) masers overlaid (note that the most central OH maser in the main image is in fact two OH sources 
close together as shown in the second sub-image). Also overlaid are the water maser source numbers as shown in Table [T] The positions 
of the methanol and OH masers have also been obtained from ATCA observations and have a similar positional accuracy to the water 
maser positions l|Caswelli 1998: EUingscn 2005). In the main image, the size of the shapes is much larger than the positional accuracy, 
however, the four sub-images (of water masers 9, 4, 3 and 2 respectively) show the positions of the maser species and indicate the 
positional accuracy of these masers overlaid on a S.O-fim GLIMPSE image of the region. Also overlaid on the sub-images are the 22-GHz 
continuum contours detected in our ATCA observation. The first contour in each case is at the 5a- level for the continuum image and 
they increase in factors of \/2. Details of the continuum sources can be found in Table |3] 
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Table 1. 22-GHz water masers detected within the survey regions. Column 1 gives the water maser source number (which is used in 
later tables), column 9 indicates whether or not each maser was detected in the ATCA observation and column 10 gives the water maser 
references. ATCA positions are quote d for all water masers w ith th e exception of G 333.29-0.38 (sourc e number 6). R eferences: * = new 
source; 1 = ICaswell et al.l 2 = lKaufmann et al.l (Il976l) : 3 = lBraz fc Scalisei ( 1982 ) 4 = Johnston et al.l lll972h . 



Source 


Water 


Right 


Declination 


Peak 


Peak Vel. 


Velocity 


epoch 


ATCA 


Ref 


number 


maser 


Ascension 


(J2000) 


Flux 


wrrt LSR 


Range 




detection? 






{l,b) 


(J2000) 




(Jy) 


(km s-i) 


(km s-i) 








1 


G 332.653-0.621 


16:19:43.569 


-51:03:37.06 


28.9 


-45.3 


-59,-43 


2006 July 


yes 


2 


2 


G 332.826-0.549 


16:20:11.089 


-50:53:16.07 


239.4 


-59.1 


-69,-45 


2006 July 


yes 


3 


3 


G 333.060-0.488 


16:20:58.002 


-50:40:46.32 


64.3 


-8.7 


-13,5 


2006 July 


yes 




4 


G 333.121-0.434 


16:20:59.762 


-50:35:51.55 


161.1 


-57.7 


-60,-46 


2006 July 


yes 


1 


5 


G 333.221-0.402 


16:21:17.913 


-50:30:17.99 


9.6 


-52.0 


-58,-48 


2006 July 


yes 




6 


G 333.29-0.38 


16:21:30.4 


-50:26:34 


7 


-49.0 


-51,-47 


2005 June 


no 




7 


G 333.364-0.358 


16:21:44.319 


-50:22:21.08 


3.2 


-52.6 


-55,-49 


2006 July 


yes 




8 


G 333.428-0.380 


16:22:07.539 


-50:20:34.29 


12.4 


4.1 


-6,6 


2006 July 


yes 




9 


G 333.608-0.215 


16:22:11.060 


-50:05:55.98 


103.3 


-49.2 


-64,-38 


2006 July 


yes 


4 



made in 1976 August by iBatchelor et al.l (|l980l ) showed a 
decline in the -51 km s~^ peak, while improved sensitivity 
uncovered additional spectral features. We observed the cur- 
rent peak intensity to be 161 Jy at -57.7 km s~^ with emis- 
sion covering the velocity range -60 to -46 km s"^ The most 
prominent feature exhibited variation by a factor of two over 
a nine-month period. 

This source is located within the RCW 106A structure 
and is offset fr om two OH masers by about 50 arcseconds 
JCaswelll Il998l). The 6.7-GHz methanol maser G 333. 121- 
0.434 ( Ellingsenll2005l ) is separated from the water maser by 
less than an arcsecond and emission is seen over a similar 
velocity range to the water maser. This maser falls within 
FIR, ^"^CO and 1.2-mm dust clumps with angular separa- 
tions of 59, 40 and 47 arcseconds from the centre of the 
respective clumps. We detected a 22-GHz continuum source 
separated from the water maser by approximately 50 arc- 
seconds whi ch may be ass ociated with the two OH masers 
observed bv lCasweUl (|l998l ) (see second sub-image of Fig. [2]). 

5. G 333.221-0.402: The peak flux density of this maser 
has remained approximately constant over the course of our 
observations, however, initial observations made in 2005 Au- 
gust showed only one spectral feature, while observations 
with the ATCA on 2006 July detected four additional peaks. 
While three of these features can be explained by the im- 
proved sensitivity offered by the ATCA the remaining sec- 
ondary peak of approximately 7 Jy should have been de- 
tected in initial observations, suggesting this maser exhibits 
some variability. 

This water maser appears to be associated with 
GLIMPSE point source G 333.2205-00.4024 and is separated 
from the centre of the nearest ^^CO clump by 33 arcseconds. 
This source falls within 10 arcseconds of the centre of a 1.2- 
mm dust clump an d 42 arcseconds from the centre of a FIR 
source detected bv lKarnik et all ((iooj). 

6. G 333.29-0.38: This water maser was discovered on 
2005 June 26 when it was detected in both polarisations and 
in adjacent spectra, however, when follow-up observations 
were made during 2005 September and October the peak 
flux density of the maser was less than 1 Jy. The earlier ob- 
servations showed the maser to have a single spectral feature 
of approximately 7 Jy at -49 km s~^, co mparable to the ve - 
locity of the Hll region, GAL 333.3-00.4, (|Huang et al.lll999l ') 



which is situated 17 arcseconds away and has a velocity of 
-52.1 km s-^ The M5X source G 333.2898-00.3898 may be 
associated with the maser as the sources are separated by 
33 arcseconds, less than the maser positional uncertainty of 
1 arcminute. At the time of the observations made with the 
ATCA the peak flux density of this maser was less than the 
5(7 detection limit of 0.75 Jy. 

7. G 333.364-0.358: A decline in the peak flux density 
has been observed since this maser was discovered on 2005 
August 25. In the initial observations a peak flux density of 
9 Jy was observed compared to the flnal observations made 
with the ATCA where a peak flux of 3.2 Jy was recorded. 
This source is separated from the centre of the nearest ^"^CO 
clump by 47 arcsecond s and t he centre of a FIR source de- 
tected bv lKarnik et al.l (|200ll ) by 37 arcseconds. 

8. G 333.428-0.380: This maser consists of two main 
spectral features at velocities of -4.5 and 4.1 km s~^ with 
flux densities of 4 Jy and 12.4 Jy respectively. A decrease 
in the flux density of the primary feature has been observed 
since observations made during 2005 April when it had a 
peak flux density of 26 Jy. This is the only maser detected 
in our survey for which there is no IRAS, MSX, 1.2-mm dust 
or other maser species within 2 arcminutes. There is however 
a nearby "CO emission peak located at about 40 arseconds 
from the maser. Like water maser G 333.060-0.488, the ve- 
locity of this maser indicates that it probably belongs to 
another molec ular cloud locate d at a distance of approxi- 
mately 1 kpc (iBains et al.ll2006l ) (see Section [3. l.ip . 

9. G 333.608-0.215: This is one of the earliest wa- 
ter masers to be di scovered and was flrst detected by 
Ijohnston et al.l (Il972l ) who observed a primary feature at 
-49 km s~^ and a secondar y feature at -57 km s~ ^. Obser- 
vations in 1976 August by IBatchelor et al] (|l980l ) found a 
100 Jy peak at -52 km s~^. We detected emission over a 26 
km s~^ velocity range with the most intense feature of 103 
Jy at -49.2 km s~^ and a decline in the -52 km s~^ feature 
to about 38 Jy. T here is an associated OH maser observed 
bv lCaswelil (|l998l ) which is displaced from the water maser 
by less than an arcsecond and shows emission over a veloc- 
ity range of -48 km s~^ to -36km s~^ (see flrst sub-image 
of Fig.lDl. 

This maser is associated with the well known Hll region 
G 333.6-0.2 which is almost totally obscured at optical wave- 
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Figure 3. Three colour GLIMPSE image of the GMC where red=8.0-/im green=5.8-/im and blue=3.6-/im. The positions of the nine 
water masers detected in this survey are represented by circles, posit ions of the me thanol masers observed by lEUingsenl ||2005|) are 
represented by squares and the positions of the OH masers observed bv lCaswelll l ll998l) are represented by crosses. Water maser sources 
1-9 are numbered in order of increasing galactic longitude. 



lengths but is one of the brightest objects at longer wave- 
lengths. This maser is probably associated with an IRAS 
source, FIR source, ^'^CO clump, CS emission and 1.2-mm 
dust clump. 

3.2 22-GHz continuum sources 

Four radio continuum sources were detected in the ATCA 
observations. Their properties are summarised in Table [3] 
and are shown in the sub-images of Fig.|2] As there is no high 
resolution, low frequency continuum data available for the 
region, spectral indices (and therefore likely source mecha- 
nisms) could not be determined. 

G 333.826-0.549 This is the strongest radio continuum 
source detected in our ATCA observations and its peak is 
located within 2 arcseconds of the water maser G 332.826- 



0.549. The nearest infrared source is /i^AS 16164-5046, with 
which it is possibly associated. 

G 333.060-0.489 This is the weakest continuum source 
that we detected. The source is spatially coincident with 
GLIMPSE point source G 333.0600-00.4888 as weU as the 
water maser G 333.060-0.488. As the separation between the 
continuum source and the near-by water maser is less than 
an arcsecond it seems likely that the two are associated. As 
this water maser (source number 3) is probably associated 
with a foreground cloud (see Section [3. l.ip and not the main 
GMC, it is likely that this continuum source also belongs to 
the foreground cloud. 

G 333. 135-0.432 This continuum source, while not asso- 
ciated with any water masers, do es appear to be associated 
with two OH masers observed bv lCaswelil (| 19981 ) (shown in 
Fig. E}. There is no associated infrared source. 
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Figure 4. Spectra of the 22-GHz water maser sources discovered in this survey. All of the spectra presented are from the ATCA 
observations with the exception of G 333.29-0.38 which is from the Mount Pleasant observations. 



Table 2. Summary of all possible associations. Details of GLIMPSE, IRAS, MSX, FIR, 1.2-mm dust, ^^CO and CS sources can be seen 
in TablesU [S] [6] 1101 [7]and |9] y indicates the presence of an association, n indicates that there is n o association , * ind icates the ^^CO 
emission peaks that we identify that were outside the velocity range of clump analysis carried out bv I Bains et al. I ll200d) and - indicates 
that no data over a similar velocity to the respective water masers was available. 
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Figure 5. ATCA spectra of the 22-GHz water maser sources detected in the search that have been previously discovered. 



Table 3. 22-GHz continuum sources. Column 1 is the source number of the nearest water maser, Columns 2-4 give the position of the 
continuum source, column 5 gives the peak of the continuum source in mjy/beam, column 6 gives the total flux density of the continuum 
source in mjy and column 7 gives the angular separation between the continuum source and the nearest water maser source. 



Source 


Continuum 


Right Ascension 


Declination 


Fpeak 


Total 


Separation 


number 


Source 


(J2000) 


(J2000) 


(mJy/beam) 


Flux Density 


from maser 




0,6) 








(mJy) 


(arcsec) 


2 


G 332.826-0.549 


16:20:11.089 


-50:53:14.07 


948 


1180 


2 


3 


G 333.060-0.489 


16:20:58.002 


-50:40:47.32 


91 


126 


1 


4 


G 333.135-0.432 


16:21:03.017 


-50:35:12.54 


720 


901 


50 


9 


G 333.605-0.212 


16:22:09.605 


-50:05:59.98 


569 


631 


15 



G 333.605-0.212 This continuum source is coincident 
with the MSX source G 333.6046-00.2124 and IRAS source 
/TJAS 16183-4958. These infrared sources are separated from 
the the peak of the continuum emission by 1.3 and 5.8 arc- 
seconds respectively. This source is offset from the nearest 
water maser that we observe by about 15 arcseconds. 



4 DISCUSSION 

4.1 Association with other maser species 

To our knowledge, this is the first high-mass star forma- 
tion complex for which untargeted searches have been made 
in all of main-line OH, 6.7-GHz methanol and 22-GHz wa- 
ter masers. Previous targeted searches for water masers 
towards known 6.7-GHz meth anol masers achieved detec- 
tion r ates of around 50 percent l|Szymczak. Pillai fc MentenI 
I2OO5I ). In contrast we find that only 25 percent of the 6.7- 



GHz methanol masers that fall within our survey region have 
an an associated water maser. We additionally find that only 
11 percent of the water masers that we detect have an as- 
sociated methanol maser which implies that targeted water 
maser searches towards 6.7-GHz methanol masers may not 
be the most efficient way to increase the number of known 
water masers. It also supports the hypothesis that water 
masers may be the most prevelant species within star forma- 
tion regions as our relatively insensitive survey has detected 
twice as many water masers as either 6.7-GHz methanol or 
OH in the corresponding regions. 



4.2 Association with infrared sources 

lEUingsenl l|2006l) found that approximately two-thirds of 6.7- 
GHz methanol masers have an associated GLIMPSE point 
source, and less than 10 percent of sources are not associ- 
ated with mid-infrared emission (at the sensitivity of the 
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GLIMPSE observations). A search of the GLIMPSE cata- 
logue reveals four of the water maser sources detected (all of 
them new discoveries) have an associated GLIMPSE point 
source within 2 arcseconds. The details of these GLIMPSE 
point sources are summarized in Table U Of the remaining 
five water masers, the four previously detected sources are all 
clearly projecte d against regions of bright mid-infrared emis- 
sion (see Fig. [3)) . ' EllingsenI (|2006^ suggested that a search for 
6.7-GHz methanol masers towards GLIMPSE point sources 
meeting the criteria [3.6]-[4.5] > 1.3 mag and 8.0 fim mag- 
nitude < 10 would detect more than 80 percent of this class 
of maser. From Table [3] it can be seen that three of the 
four GLIMPSE point sources associated with water masers 
satisfy these criteria. We compared the [3.6]-[4.5] colours of 
the GLIMPSE sources associated with water masers with 
the colours of GLIMPSE sources associated with 6.7-GHz 
methanol masers (see fig. 16 of Ellingscn). The colours of 
the sources associated with water masers are clustered in 
the less red end of the range observed in the methanol 
associated sources. The idea that different maser species 
may trace different phases of the high-mass evolutionary se- 
quence is not new, however, it has been receiving renewed 
attention lately, as sensitive, high-resolution observations 
at sub-millimetre thro ugh mid-infrared wavelengths become 
more readily available. lEllingsenI l|2006l ) looked at this ques- 
tion in some detail and we will not repeat the arguments 
here. However, the striking difference between the relative 
location of the water and methanol masers in the G 333.2- 
0.6 GMC, and the less red colours of the water maser as- 
sociated GLIMPSE sources supports the hypothesis that 
6.7-GHz methanol masers may trace a generally earlier evo- 
lutionary phase than water masers. The very small water 
maser sample size prevents us from drawing any firm conclu- 
sions, however, it suggests that comparison of the properties 
of GLIMPSE sources associated with water and methanol 
masers may provide useful insights relating to this question. 

The top and bottom inset images in Fig. [2] each show a 
ring of diffuse 8.0-fim emission surrounding a darker region. 
We have closely examined the GLIMPSE images of these 
regions. In the case of the southern most region (associated 
with the water maser G 333.826-0.549) it is clear that this 
is the result of a known image artifact that is present in 
some GLIMPSE sources where saturation occurs. An expla- 
nation for the northern most source (associated with the 
water maser G 333.608-0.215) is not as clear. The ring like 
structure seen in the 8.0-/im is similarly seen in the other 
GLIMPSE bands with radius decreasing with wavelength 
and the point where the dip in fiux density per steradian 
occurs is close to the GLIMPSE saturation level. Interest- 
ingly, the MSX images of the region show a similar structure, 
warranting further observations of this source at both radio 
and infrared wavelengths to better understand its nature. 

Three of the nine water masers detected in this survey 
(all previously known sources) have an IRAS source within 
30 arcseconds. Two of the water masers that we observe have 
an associated MSX source. The details of these IRAS and 
MSX sources are presented in Table [5] 

iKarnik et al.l (|200ll ) detected 23 FIR sources in a survey 
of the GMC at 150 and 210 /xm. Of the six most luminous of 
these sources, five have associated water masers (see Table[6] 
for details). 




-120 -100 -80 -60 -40 -20 20 



Velocity [km/s] 

Figure 6. The mean velocity profile of the ^ '^CO emission aver - 
aged over the entire field of observations from lBains et ahl ll2006h . 
Overlaid are the velocity ranges of the water masers. (1) indicates 
the total velocity range of the seven water masers with compara- 
ble velocity ranges; maser sources 1, 2, 4, 5, 6, 7 and 9, while (2) 
indicates the velocity ranges of maser sources 3 and 8. 



4.3 Association with molecular gas and dust 
clumps 

4-3.1 CO clump associations 

I Bains et all l|2006t l produced clumpfind fits of the GMC 
using the 2D '^^CO data integrated over the velocity range 
-65 km s"'^ to -35 km s~^ which covers the most prominent 
feature of the velocity profile seen in Fig.[6l Seven of the nine 
water maser s that w e detect fall within the velocity range 
analysed by iBains et al. (2006) and all of these masers fall 

within a ^^CO clump that they identi^^ 

We obtained the "CO data of iBains etahl (|2006l ) in 
the velocity range -120 km s~^ to 20 km s~^ which has 
allowed us to identify additional ^'^CO emission peaks (or 
possible clumps) which may be associated w ith the two 
masers not coincident with a clump identified bv lBains et al.l 
(l2006l ). G 333.060-0.488 and G 333.428-0.380. These emis- 
sion peaks, like their associated water masers, have veloc- 
ity ranges which fall outside that analysed by iBai ns et al.l 
(2006). Water masers G 333.060-0.488 and G 333.428-0.380 
show emission over a comparable velocity range of a sec- 
ondary feature seen in the ^^CO emission (Fig. 

Details of the ^^CO clumps identified by iBains et al.l 

(j2006l ) and the "CO emission peaks that we identify along 
with their associated water maser sources are summarized 
in Tab le [71 while clump properties identified bv lBains et al.l 
(! 2006h are presented in Table [S] The velocity ranges of the 
^■^CO clumps or emission peaks have been obtained from 
visual inspection of the data. The '^^CO emission peaks 
have more extensive velocity ranges than the associated wa- 
ter maser emission in all cases with the exception of water 
masers G 333.060-0.488 and G 333.428-0.380. Interestingly 
these two water masers have similar velocity ranges which 
place them both at a distance of around 1 kpc (see Sec- 
tion (STTTlJ. In the majority of cases the velocity of the peak 
fiux density of the water masers are within a few km s~^ of 
the velocity of the associated peak ^^CO emission. 
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Table 4. Water masers with GLIMPSE point source associations, here IRAC bands 1, 2, 3 and 4 correspond to 3.6-/im, 4.5-^m, 5.8-^m 
and 8.0-^m. Column 1 gives the water maser source number (see Table [T1 for details), column 2 is the associated GLIMPSE point source, 
column 3 gives the angular separation between the GLIMPSE point source and the water maser source and colums 4 to 7 show the 
magnitudes of IRAC bands 1, 2, 3 and 4 for each of the GLIMPSE point sources, while column 8 gives the [3.6]-[4.5] colour for each. 



Source GLIMPSE Separation Magnitude (mag) [3.6]-[4.5] 

number point source (arcsec) IRAC IRAC IRAC IRAC colour 

{I, b) band 1 band 2 band 3 band 4 (mag) 



3 


G 333.0600-00.4888 


1.3 


11.964 


10.633 


9.258 




1.331 


5 


G 333.2205-00.4024 


1.8 


8.973 


7.434 


6.368 


5.895 


1.539 


6 


G 333.3639-00.3574 


0.8 


11.263 


9.780 


8.560 


7.551 


1.483 


7 


G 333.4285-00.3809 


2.0 


14.233 


13.373 






0.860 



Table 5. Possible infrared and water maser source associations. Column 1 gives the water maser source number (see Table[T]for details), 
column 2 shows the closest IRAS source within 2 arcminutes of the water masers, column 3 gives the angular separation between the 
water maser sources and the IRAS source, column 4 is the closest MSX source to within 1 arcminute and column 5 gives the angular 
separation between the MSX source and water maser source. 



Source 


IRAS 


Separation 


MSX 


Separation 


number 


source 


(arcsec) 


source 


(arcsec) 


1 


16158-5055 


11 






2 


16164-5046 


29 


G332.8269-00.5489 


5 


6 






G333.2898-00.3898 


33 


9 


16183-4958 


20 







4-3.2 Associations with CS emission 

The GMC has been mapped in CS emission by the Delta 
Quadrant Survey Team as part of the ongoing project at 
University of New South Wales. CS emission is sensitive to 
higher densities and as a result the mapped region is less 
extensive than that covered in the ^'^CO mapping (and has 
been obtained in advance of publication as a result of cor- 
respondence with Nadia Lo from UNSW). This means that 
CS data was not available for two of the positions of the de- 
tected water masers. We inspected the data over a velocity 
range of -80 km s~^ to -20 km s~^ and identified CS emis- 
sion peaks near all of the seven water masers for which CS 
data was available. The details of these emission peaks and 
their associated water masers are summarized in Table [9] 

Unlike the ^=*C0 emission, CS emission is observed over 
a smaller velocity range than the associated water maser 
emission, however like the ^^CO emission the velocity of the 
peak flux density of the water masers correlates to within a 
few km s~^ with the velocity of the CS emission peaks. The 
association of the water masers with CS emission indicates 
that the water masers are probably associated w ith massive 
star formation l|Bronfman. Nvman fc Mavll 19961 ) . 

4-3.3 Association with 1.2-mm dust clumps 

Seven of the water masers that we detect ar e associated with 
a 1.2- mm dust emission peak observed by iMookeriea et al.l 
l|2004l) and the propert ies of th e dust clumps are summarized 
in Table [10] Mookori ea et al.l (|2004) identifled 95 1.2-mm 
dust clumps within the GMC and 73 of these fall within 
our observed regions. A statistical analysis of the 1.2-mm 
clumps that fall within the observed regions is presented in 
Section l44l 



4.4 CO and 1.2-mm dust clump analysis 

In order to investigate the properties of the molecular gas 
and dust in the regions with associated water maser emission 
wc have fitted a Binomial generalized linear model (GLM) 
(McCuUagh & Ncldcr 1989) to the maser presence/absence 
da ta using ^"^CO and 1 .2-m m dust clump properties reported 
by iBains et al.l (|2006D and IMookeriea et al.l l|2004l ) respec- 
tively, as predictors. In the case of the ^^CO clumps the 
properties considered were the integrated flux densities of 
the clump peaks (10 K km s~^), the peak flux in terms 
of antenna temperature (K), clump radius (pc), ^''CO col- 
umn density (10^^ cm~'^) and the total LTE molecular mass 
calculated from the ^"^CO data. In the case of the 1.2-mm 
dust clump analysis the properties Fp^ak (mjy/beam), ra- 
dius (pc), total integrated flux density (F„) in Jy, mass (M0) 
and nH2 (lO'* cm~'^) where investigated as predictors. A Bi- 
nomial GLM predicts the probability, pi , of finding a maser 
in the i*'' clump, in terms of the clump properties xu X2i 
X3i . - - Xmi . The model takes the form 

yi ~ Bin(l,pi) 

Pi 

log- = Po + PlXli + f32X2i + - - - + l^mXmi 

I- Pi 

where yi is the maser presence or absence in the i"" clump 
and /3o , Pi, 132 - - - Pm are the regression coefficients to be 
estimated. 

To test the significance of individual clump properties 
all possible single term models were fitted, and compared by 
analysis of deviance to the null model consisting of only an 
intercept. Stepwise model selection based on the Akaik e In- 
formation Criteria (AIC) (jBurnham fc Andersonl[2002h was 
used to select the most parsimonious model with the greatest 
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Table 6. Separation from FIR sources, detection at 150 and 210 /^m. Positions given are of the p eak of the 210 t im |Karnik et aI.ll2001^ . 
Columns are as follows; 1 is the maser source number, 2 is the FIR source number quoted by iKarnik et al. and 4 are the 

Right Ascension and Declination of the clumps, 5 gives the angular separation between the water maser and the peak of the FIR source, 
columns 6 and 7 give the FIR source's peak flux densities at 150 and 210 fim respectively and column 8 gives the source luminosity. 



Source 


FIR 


Right 


Declination 


Separation 


150 fj.m 


210 /tm 


Luminosity 


number 


source 


Ascension 


(J2000) 


(arcsec) 


Flux Density 


Flux Density 


(10^ Lq) 




number 


(J2000) 






(Jy) 


(Jy) 




2 


Sll 


16:20:06.0 


-50:53:19 


48 


15380 


7790 


411 


4 


S15 


16:20:57.9 


-50:34:55 


59 


20811 


10009 


482 


5 


S18 


16:21:20.3 


-50:29:43 


42 


2180 


1043 


115 


6 


S20 


16:21:29.8 


-50:25:57 


37 


15572 


8015 


460 


9 


S23 


16:22:08.9 


-50:05:28 


35 


27841 


11555 


921 



Table 7. Water maser sources with nearby ^^CO clumps identified bv lBain s et al.l | |2006^ ( over velocity range -65 km s^^ to -35 km s~^ 
) or emission peaks as identified by visual inspection of the data used bv lBains et al.l l |2006l V Column 1 is the water maser source number 
(see Table [T] for details), colum n 2 and 3 give the velocity range and the velocity of the peak emission of the water maser, column 4 gives 
the '^^CO clump number from lBains et aP l|2006l ) if applicable, columns 5 and 6 show the right asc ension and declination of the ^^CO 
clump or alternatively the ^^CO emission we identified from the ^'^CO data from lBains et al JHooi), columns 7 and 8 give the velocity 
rang e and the ve l ocity of the ^'^CO emission peak as we have identified and column 9 is the separation between the centre of the clump 
that iBains et al. 1 1I2OO6I ') reported or the emission peak that we identified. The velocity ranges of the additional emission peaks that we 
identify as being asso ciated with water masers 3 and 8, indicate that they are probably located along the line of sight at a distance of 
approximately 1 kpc l lBains et al.ll200^) . 



Source 


Velocity 


Peak 




Right 


Declination 


Velocity 


Peak 


Separation 


number 


Range 


Velocity 


source 


Ascension 


(J2000) 


Range 


Velocity 


(arcsec) 




(km s-i) 


(km s-i) 


number 


(J2000) 




(kms-i) 


(kms^^) 




1 


-59,-43 


-45.3 


9 


16:19:40 


-51:03:27 


-60,-43 


-50 


35 


2 


-69,-45 


-59.1 


7 


16:20:11 


-50:53:27 


-65,-44 


-55 


11 


3 


-13,5 


-8.7 




16:20:58 


-50:40:39 


-9,-2 


-6 


8 


4 


-60,-46 


-57.7 


1 


16:21:03 


-50:35:27 


-64,-45 


-51 


40 


5 


-58,-48 


-52.0 


14 


16:21:21 


-50:30:03 


-57,-40 


-52 


33 


6 


-58,-47 


-49.0 


2 


16:21:31 


-50:26:51 


-58,-45 


-52 


23 


7 


-55,-49 


-52.6 


22 


16:21:42 


-50:21:39 


-67,-42 


-51 


48 


8 


-6,6 


4.1 




16:22:05 


-50:21:02 


-1,10 


6 


40 


9 


-64,-38 


-49.2 


5 


16:22:08 


-50:06:27 


-61,-32 


-46 


42 



predictive power. The AIC is a trade off between goodness 
of fit and model complexity and is defined as 



AIC — — 2(max log likelihood) -|- 2(number of parameters) 



with the preferred model bein g the model with the lowest 
AIC jVenables fc Riplevll2002l ). 

For ease of comparison between the data sets box plots 
were created for each of the clump properties. The solid hor- 
izontal line in each of these plots represents the median of 
the data. The box represents the 25*'' to the 75*'' percentile, 
while the vertical line from the top of the box represents 
data from the 75*'' percentile to the maximum value and 
the vertical line from the bottom of the box represents data 
from the 25*'' percentile to the minimum value. Outliers are 
represented by dots. Box plots of each of the clump proper- 
ties of both the ^^CO clumps and the 1.2-mm dust clumps 
are shown in Figs. [7] and |8] respectively. 



CO clump results 

Fits of the single term additi on Binomial model to the ^"^CO 
clump properties reported bvlB ains et al. 1(12006'), showed an 
increasing probability of the presence of a maser was associ- 
ated with all of the tested factors (the integrated flux density 
of the clump peak, ^^CO clump peak in terms of antenna 
temperature, clump radius, ^^CO column density and the 
total LTE molecular mass calculated from the ^"^CO data). 
This means that any of the clump properties (in isolation) 
gives an indication of the likelihood of maser presence. Ta- 
ble [TT] gives a summary of the single term addition Bino- 
mial model. The same information is shown graphically in 
Fig. [7] which clearly illustrates that for all clump properties 
there is a difi'erence between the ^^CO clumps that have an 
associated water maser and those that do not. In general 
the ^"^CO clumps with associated water masers are bigger, 
denser, brighter and more massive. The most parsimonious 
model for predicting maser presence involved the integrated 
flux density of the clump peak, clump radius and the total 
LTE molecular mass calculated from the ^^CO data. This 
means that if the integrated flux density of the peak, radius 
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Table 8. Properties of the j^CO clumps llBaina et al.ll2006h . Column 1 is the number of the associated water maser source, column 2 
is the ^"^CO clump number l l Bains et al.ll2006l) r column 3 is the integated flux density of the clump peak, column 4 is the ^'^CO clump 
peak in terms of antenna temperature, column 5 is the clump radius, column 6 is the ^'^CO column density and column 7 is the total 
LTE molecular mass. 



Source ^-^CO Integrated Flux Peak Radius N(i^CO) Mass 
number Clump (10 K km s^^) (K) (pc) (IQi"^ cm^^) (lo^ Mq) 



1 


9 


10.5 


13.3 


1.7 


6.6 


4.6 


2 


7 


12.0 


14.4 


1.8 


6.2 


5.0 


4 


1 


18.5 


46.0 


2.6 


10.1 


16.0 


5 


14 


8.8 


15.7 


1.8 


6.9 


5.5 


6 


2 


15.8 


41.4 


2.4 


10.3 


14.4 


7 


22 


8.9 


25.0 


2.2 


7.6 


8.7 


9 


5 


13.8 


27.4 


2.2 


7.9 


9.5 



Table 9. Water masers sources with nearby CS emission peaks. Column 1 is the water maser source number, columns 2 and 3 give the 
velocity range and velocity of emission peak, columns 4 and 5 are the right ascension and declination of the CS emission peak, columns 
6 and 7 show the velocity range and the emission peak of the CS emission and column 8 shows the angular separation between the water 
maser source and the CS emission peak. 



Source 


Velocity 


Velocity 


CS 


Clump 


Velocity 


Peak 


Separation 


number 


Range 


Peak 


Right Ascension 


Declination 


Range 


Velocity 


(arcsec) 




(kms-i) 


(kms-i) 


(J2000) 


(J2000) 


(kms-i) 


(kms-i) 




1 


-59,-43 


-45.3 


16:19:40 


-51:03:30 


-56,-44 


-50 


33 


2 


-69,-45 


-59.1 


16:20:11 


-50:53:19 


-61,-51 


-55 


4 


4 


-60,-46 


-57.7 


16:21:03 


-50:34:56 


-61,-44 


-53 


63 


5 


-58,-48 


-52.0 


16:21:22 


-50:30:46 


-57,-48 


-51 


45 


6 


-58,-47 


-49.0 


16:21:29 


-50:26:31 


-58,-45 


-51 


12 


7 


-55,-49 


-52.6 


16:21:41 


-50:23:19 


-55,-45 


-49 


69 


9 


-64,-38 


-49.2 


16:22:08 


-50:06:17 


-58,-39 


-47 


40 



Table 1 1 . Analysis of deviance table for the single term models 
(using ^^CO clump properties), showing the AIC and the de- 
viance together with the associated likelihood ratio statistic and 
p-value for the test of the hypothesis that the stated single term 
model provides no better fit than the null model consisting only 
of an intercept. 



Predictor 


AIC 


Deviance 


LRT 


p-value 


none 


39.098 


37.098 






Integrated 


28.774 


24.774 


12.325 


0.000447 


Peak 


33.093 


29.093 


8.005 


0.004664 


Radius 


34.063 


30.063 


7.035 


0.007994 


Density 


34.615 


30.615 


6.483 


0.010892 


Mass 


35.265 


31.265 


5.833 


0.015729 



and mass is known for a CO clump then a probability of 
maser presence can be determined. The estimated regression 
relation is 

log— ^ = -21.2018 + l.2,Qi^X^r^tegrate.d + 8.0589x,ad,„, 
I- Pi 

1 .2290x mass 5 

where Xintegrated, Xradius and Xmass represent clump proper- 
ties integrated flux density of the clump peak (10 K km s~^), 
radius (pc) and mass (10"^ Mq). The full regression summary 
is shown in Table [12] 



Table 12. Summary table for the Binomial regression model, 
showing for each predictor the estimated coefficient and its stan- 
dard error, and the standardised z-value and p-value for the test 
of the hypothesis that /3i=0. 



Predictor 


Estimate 


Std. Error 


z value 


p-value 


Intercept 


-21.2018 


9.0538 


-2.342 


0.0192 


Integrated 


1.3037 


0.5041 


2.586 


0.0097 


Radius 


8.0589 


5.0069 


1.610 


0.1075 


Mass 


-1.2290 


0.6153 


-1.997 


0.0458 



Setting the threshold probability of the ^^CO clump 
model at 0.5 (i.e. a value greater than 0.5 suggests a water 
meiser will be associated with a clump, while a lower value 
suggests no water maser) , we flnd the misclassification rates 
to be low. Of the 40 ^''CO clumps that fall in our observa- 
tion regions which correspond the the main part of the GMC 
(—65 km s~^< VLSR < —35 km s~ ), seven have an associ- 
ated water maser while our model predicts that five of these 
have an associated water maser and returns a false-negative 
for the remaining two clumps. There are 33 clumps within 
the survey regions that do not have associated water masers, 
our model predicts that 31 of these do not have associated 
water masers but returns a false-positive for the remaining 
two clumps. 
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Table 10. Maser sources with nearby 1.2-mm dust clumps identified by iMookeriea et al ] l|2004l '). Column 1 is the water maser source 
number, column 2 is the dust clump number, columns 3 and 4 give the right ascension and declination of the dust clump, column 5 is 
the angular separation between the water maser and associated dust clump, column 6 gives the peak flux density for the clump, column 
7 gives the clump radii, column 8 is the total integrated flux densities of the clumps, column 9 gives the estimated mass of the clumps 
assuming a temperature of 40 K and column 10 gives the column density. 



Source 
number 


1.2-mm 

dust 
source 


Right 
Ascension 
(J2000) 


Declination 
(J2000) 


Separation 
(arcsec) 


F , 

peak 

(mJy/beam) 


Radius 
(pc) 


(Jy) 


Mass 
(Mq) 


^H2 

(10" cm-3) 


1 


MMS84 


16:19:38.9 


-51:03:28 


45 


2499 


1.48 


21.9 


2871 


2.45 


2 


MMS68 


16:20:11.9 


-50:53:17 


8 


12460 


1.70 


36.1 


5548 


2.56 


3 


MMS51 


16:20:52.1 


-50:40:51 


56 


348 


1.11 


2.7 


1029 


0.73 


4 


MMS39 


16:21:03.7 


-50:35:23 


47 


8925 


1.19 


32.6 


4584 


6.18 


5 


MMS33 


16:21:18.6 


-50:30:25 


10 


1095 


0.90 


3.5 


1325 


1.76 


6 


MMS29 


16:21:32.7 


-50:27:12 


44 


6502 


1.27 


24.7 


9281 


4.38 


9 


MMS5 


16:22:10.1 


-50:06:06 


14 


40013 


1.88 


129.5 


15936 


5.44 



Table 13. Analysis of deviance table for the single term models 
(using the 1.2-mm dust clump properties), showing the AIC and 
the deviance together with the associated likelihood ratio statistic 
and p- value for the test of the hypothesis that the stated single 
model provides no better fit than the null model consisting only 
of an intercept. 



Predictor 


AIC 


Deviance 


LRT 


p- value 


none 


48.130 


46.130 








Fpeafc 


28.088 


24.088 


22.042 


2.668e- 


-06 


Radius 


24.330 


20.330 


25.800 


3.787e- 


-07 




25.428 


21.428 


24.702 


6.693e- 


-07 


Mass 


27.467 


23.467 


22.663 


1.930e- 


-06 


Density 


39.969 


35.969 


10.160 


0.001435 



l-2-mm dust clump results 

In terms of the 1.2-mm dust clump properties, fits of the 
single term addition Binomial model showed an increas- 
ing probability that the presence of a maser was associated 
with increasi ng valu e of all of the clump properties reported 
by rMookc rica et al.l (|2004). Table [13] gives a summary of 
the single term addition. This means that any one of the 
clump properties may give an indication of the likelihood of 
maser presence. There is a significant difference between the 
clumps with associated water masers and those without, evi- 
dent from Fig.[S] Clumps with associated masers are bigger, 
denser, more massive and have higher flux densities than 
clumps where we see no associated water maser. The most 
parsimonious model for predicting maser presence involved 
only the radius of the 1.2-mm dust clumps. This equation al- 
lows the probability of maser presence to be predicted know- 
ing only the radius of a 1.2-mm dust clump. The estimated 
regression relation is 



log-; 



= -11.477 -f9.174a;^ 



where Xradius is the radius of the 1.2-mm dust clump in pc. 
The regression summary of this model is shown in Table [T4l 
The misclassification rates for the model, given a proba- 
bility of 0.5 of maser presence within a given clump are again 
low in predicting the clumps that have no associated water 



Table 14. Summary table for the Binomial regression model, 
showing for each predictor the estimated coefficient and the stan- 
dardised z-value and p-value for the test of the hypothesis that 
ft=0. 



Predictor 


Estimate 


Std. Error 


z value 


p-value 


Intercept 


-11.477 


3.537 


-3.245 


0.00118 


Radius 


9.174 


3.163 


2.900 


0.00373 



maser. Of the 73 1.2-mm dust clumps that fall within our 
survey region, seven have an associated water maser, while 
66 do not. Our model predicts that 65 of the 66 clumps with 
no associated water maser will not have an associated water 
maser and returns a false-positive result for the remaining 
clump. The model predicts that of the seven clumps that 
have an associated water maser only four will have an as- 
sociated water maser and returns a false-negative for the 
remaining three clumps. 

As this model is only concerned with clump radius it 
is easy to illustrate its physical implications. By setting the 
probability of maser presence to be 0.5, for example, we 
are able to determine that the corresponding clump radius 
is approximately 1.25 pc. This means that all 1.2-mm dust 
clumps with a radius of 1.25 pc or greater have a probability 
of 0.5 or more of having an associated water maser. 



5 CONCLUSIONS 

Regions within the GMC associated with RCW 106 have 
been surveyed for 22-GHz water masers. This resulted in 
the detection of nine water masers (five of these being new 
detections) and four 22-GHz continuum sources. All of the 
water masers that we observed have exhibited some level 
of variability over the 11 month course of these observa- 
tions. The most extensive temporal variability was observed 
in water maser G 333.060-0.488, which showed a variation 
in peak flux density of a factor of 10 over the observa- 
tional period. The GMC has previously been s earched for 
6.7-GHz methano l masers llEllingsen et al" Il996l) and main- 
line OH masers l|Caswell. Havnes &: Goss 1980l ). In addi- 
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Figure 7. Box plots of each of the ^^CO clump properties split 
into the categories of yes and no, referring to maser presence and 
absence respectively. 



Figure 8. Box plots of the 1.2-mm dust clump properties split 
into the categories of yes and no, referring to maser presence and 
absence respectively. 



A search for 22- GHz water masers 17 



tion to the nine water masers detected ther e are four 6.7- 
GHz methanol masers jElhngsen et al. | ll996l ) and three OH 
masers (jCaswell. Havnes fc Gosd ll98C ) within the regions 
surveyed here. AU of the three species of masers have sub- 
arcsecond positional accuracy which allows a comparison of 
the relative positions of the respective maser species. The 
water masers that we detect lie along the main axis of star 
formation within the GMC while the methanol masers are 
located near the periphery. We find there to be little overlap 
between the sites of the different maser species, in fact there 
are only two associations with other maser species. There is 
one association with a 6.7-GHz maser and one association 
with an OH maser. 

Four of the new water maser detections are associated 
with GLIMPSE point sources, of similar colours to those 
associated with detected 6.7-GHz methanol masers. There 
is a slight bias for the water maser associated sources to be 
less red. This coupled with the relative positions of the water 
masers and the 6.7-GHz methanol masers with respect to the 
main axis of star formation lends support to the hypothesis 
that 6.7-GHz methanol masers trace an earlier evolutionary 
phase than water masers. 

All of the water masers are associated with a ^"^CO 
emission pe a k tha t we identify or a clump reported by 
iBains et all (120061 '). Statistical investig ation of the " CO 
l|Bains et aLll2006i ) and 1.2-mm dust ( Mookeriea et al.ll2004l ') 
clumps shows that there is a strong increase in likelihood of 
water maser detection with increased clump radius, mass, 
density and brightness. After fitting a Binomial generalized 
linear model to the maser pre senc e data using the clump 
properties of lBains et"al] (|2006l ') and lMookeriea et al.1 (|2004l ') 
as predictors we obtained the simplest models with the 
greatest maser predictive power. In the case of the 1.2-mm 
dust clumps our model uses only clump radius to predict the 
likelihood of the clump having an associated water maser. 
However the model generated for ^^CO clumps takes into ac- 
count radius, integrated flux density of the clump peak and 
mass. These models have a low misclassification rate and 
may allow more efficient targeted searches for water masers 
(where appropriate ^'^CO or 1.2-mm data is available) than 
those previously conducted towards other species of masers 
or mid-infrared sources. While our survey was of compara- 
tively small scale, we believe that the results are indicative 
of the likelihood of occurrence of water masers with respect 
to "CO and 1.2 -mm dust clumps. 

Observations of ^"^CO clumps and 1.2-mm dust clumps 
within the G 333.2-0.6 GMC that do not fall within the re- 
gions already surveyed will allow us to test and refine our 
predictive models. Additional observations towards ^'^CO 
clumps and 1.2-mm dust clumps which are not part of the 
G 333.2-0.6 GMC will allow us to determine models with 
more highly accurate predictive properties and lower stan- 
dard errors. As more data for the G 333.2-0.6 GMC in other 
molecular line transitions becomes available, analysis similar 
to what we have undertaken here will be carried out. The 
multitude of data available for this GMC provides a unique 
oppourtunity which may provide the evidence to support an 
evolutionary argument. 
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